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Abstract 

We have observed nanopattern formation with robust and controllable spatial ordering by 
laser-induced dewetting in nanoscopic metal films. Pattern evolution in Co film of thickness 
1 < h < 8 ram on SiC>2 was achieved under multiple pulse irradiation using a 9 ns pulse 
laser. Dewetting leads to the formation of cellular patterns which evolve into polygons that 
eventually break up into nanoparticles with monomodal size distribution and short range or- 
dering in nearest-neighbour spacing R. Spatial ordering was attributed to a hydrodynamic thin 
film instability and resulted in a predictable variation of R and particle diameter D with h. 
The length scales R and D were found to be independent of the laser energy. These results 
suggest that spatially ordered metal nanoparticles can be robustly assembled by laser-induced 
dewetting. 
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1 Introduction 



Ordered metal nanoparticle arrays have become an active area of research due to their potential 
applications in nonlinear optics and nanophotonics below the diffraction limit [1-3]. However, 
practical realization of such applications require that cost-effective and reliable processes be de- 
veloped for the manufacturing of ordered metal nanoarrays. Physical phenomena that lead to the 
formation of patterns with characteristic length scales and predictable time scales could be ex- 
ploited in designing such nanomanufacturing processes. For instance, pattern formation following 
ion irradiation has been known to modify surface topography [4-8]. Rippling of surfaces under 
ion irradiation appears to be primarily due to an instability resulting from the competition between 
ion erosion and smoothening due to surface diffusion with the length and time scales determined 
by ion flux, surface temperature and surface diffusion parameters [9-12]. Another route to creat- 
ing nanopatterns could be through instabilities of a thin fluid film leading to spatial patterns with 
length and time scales that depend on the thermophysical material properties such as interfacial ten- 
sion, contact angle with the substrate, fluid viscosity and, for ultrathin films, long range dispersion 
forces such as the van der Waal's interaction. An example especially pertinent to thin film pat- 
tern formation is the hydrodynamic dewetting instability such as spinodal dewetting which occurs 
when attractive intermolecular forces exceed the stabilizing effect of interfacial tension [13-16]. 
Under such conditions, spontaneous film thickness fluctuations could be amplified and this will 
result in the breakup of the film leading eventually to the formation of drops/particles with well 
defined spatial order [17]. Dewetting dynamics leading to particles have been studied in detail in 
polymer films that are in the liquid state close to room temperature [18-21]. However, ordered 
nanoparticle formation via dewetting in metal films has remained relatively unexplored , primarily 
due to the high temperatures required to melt the metal and observe dewetting within practical 
time scales. Earlier work has shown that ion-irradiation of ultrathin Pt films on Si0 2 substrates 
results in dewetting patterns with spatial order [22,23]. More recently our group has shown that 
multiple ns pulsed laser melting of nanoscopic Co films can lead to short- and long-range spatial 
order [24-26] . However, the the relationship of particle spacing R and diameter D with film thick- 
ness h was not determined. Furthermore, it is not clear whether thin film hydrodynamic models 
could help predict such relationships at nanometer length scales. Another unanswered question is 
whether dewetting under multiple cycles of phase change, i.e. melting and resolidification, would 
lead to robust final states. In this article, we address the above questions quantitatively by ana- 
lyzing the results of a series of systematic experiments and interpreting them with the aid of thin 
film hydrodynamic theories. Such comprehensive understanding could help in knowledge-based 
design of nanomanufacturing facilities for ordered metal nanoarrays. 

Specifically, we investigate the evolution of dewetting patterns in nanoscopic Co metal films 
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on Si0 2 in vacuum under multiple pulse irradiation with a laser of pulse length t p = 9 ns and 
wavelength 266 nm. Co nanoparticles are useful towards nanoscale magnetic applications and in 
catalysis. Dewetting was studied as a function of film thickness h, laser energy density E and 
laser irradiation time expressed in terms of the number of pulses n. We discovered that pattern 
formation evolved from cellular structures composed of polygons at early stages, which coalesced 
and eventually changed into nanoparticles at later times. Every observed pattern was ordered with 
a characteristic length scale while the nanoparticles had a well-defined nearest-neighbour (NN) 
spacing R and a monomodal size distribution with an average diameter D. The observed variations 
of R and D as function of h were in agreement with the linear theory of spinodal dewetting. We 
also found that R was independent of n and E, provided net Co evaporation was small. Based 
on these observations we infer that once the pulsed laser melting initiated dewetting, subsequent 
pulses evolved this instability while maintaining its spatial characteristics, thereby leading to a pre- 
dictable final state. This offered strong evidence that the dewetting dynamics is extremely robust. 
To our knowledge this is the first experimental evidence for: (i) formation of metal nanoparticles 
with spatial order via a hydrodynamic dewetting instability; ( ii) evidence for the linear spinodal 
theory under multiple instances of phase change; and (Hi) evidence for robustness of the dewetting 
process to experimental conditions. 



2 Experimental Details 

Co films with thickness ~ 1 to 8 nm were deposited at rates varying between 0.5 to 2 nm/min 
onto commercial optical quality Si0 2 /Si wafers consisting of 400 nm thick thermally grown ox- 
ide layer on polished Si(100) wafers under high vacuum (2 x 10~ 8 Torr) by e-beam evaporation 
at room temperature [24, 26]. For every film thickness we measured the surface roughness and 
average grain size. Atomic force microscopy (AFM) of the surface roughness established an up- 
per limit of 0.1 ± .03 nm average rms roughness for the entire thickness range while transmission 
electron microscopy gave a average grain size of ~ 20 nm. No systematic change in roughness or 
particle size was detected with increasing film thickness. Further, no spatially ordered defects or 
other heterogeneities were observed on the native substrate surface. For the film thickness range 
investigated the laser absorption at 266 nm wavelength was uniform within the Co film (absorption 
depth ~ 11 nm). Therefore vertical temperature gradients could be neglected and the entire film 
temperature may be modeled as a uniform, time varying function that depends on the laser power 
(source term) and conduction of heat to the substrate (sink term) [27]. The film thickness was 
measured by concentration measurements of Co via energy dispersive X-ray spectrometry (EDS) 
in a scanning electron microscope (SEM). The concentration was converted into an equivalent 
thickness value by using calibration based on step-height and Rutherford backscattering measure- 
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ments of the film thickness. The films were irradiated immediately after deposition in vacuum at 
normal incidence by an unfocused laser beam of area 3x3 mm 2 at a repetition rate of 50 Hz. 
Under these conditions, the film surface was never exposed to air. Hence no oxide layer was ex- 
pected to form on the film surface. Pattern formation was studied for energy densities E between 
40 < E < 150 m J /cm 2 and for n between 10 and 10,500 pulses. For each thickness the range 
of E was chosen to be above the melt threshold. This was determined by a visible roughening of 
the metal film surface, as detected under high-resolution SEM within the longest time scale of the 
experiment (i.e. after 10, 500 laser pulses) [28]. Analytical estimates based on a 1-D heat transport 
model as well as finite element simulations showed that the lifetime r m of the Co liquid ranged 
from ~ 2 < r m < 12 ns in the entire h range studied [27]. The same model predicted cooling 
rates of the order of 10 10 K/s and a total heating plus cooling time per pulse of ~ 100 ns, which 
was much smaller than the spacing between pulses of 20 ms. Hence, during each pulse, the film 
melts and in the molten state, dewetting process occurs. Upon the cessation of the pulse, the film 
solidifies within a time scale much smaller than the interval between subsequent pulses. The max- 
imum temperature of the molten liquid above the Co melting point T m = 1768 K was estimated to 
be ~ 500 K. Under these conditions, the maximum evaporation observed after the longest n at the 
highest E was < 15%. 

3 Results 

Fig. dshows the dewetting morphology observed for irradiation at E = 120 mJ /cm 2 of a 3.8 nm 
thick Co film for increasing number of laser pulses n. This evolution sequence is typical of the 
ordered morphology observed in the thickness range of 1 < h < 8 nm. After a few pulses, the 
dewetting morphology is characterized by a cellular network of polygons (Fig. UK)- As n is in- 
creased the metal/substrate contact lines recede to the edge of the holes resulting in a network of 
coalescing polygonal features (Fig. [1)3). Continued irradiation results in the formation of nanopar- 
ticles primarily at the junctions of the polygons, as evident from the location of the particles in 
Fig. [l)B and in C and D. At every observed stage, a characteristic length scale L was present, as 
evidenced by the annular form of the power spectrum of the spatial correlations in the intensity 
variation within each pattern. The power spectrum is shown alongside the corresponding images 
in Fig. OK to D. For patterns consisting of polygons, L represents the average distance between the 
center to the vertices, while in the case of the nanoparticles it represents the average NN spacing 
R. 

A nonlinear peak fitting algorithm (Levenberg-Marquardt) was used to locate Gaussian shaped 
peaks in the radial distribution function g(k) versus the wavenumber k = 1/L for each pattern 
(Fig. |2j\). The variation of L with n for various E values is shown in Fig. |2jB (the lines are only 
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guides for the eye). The initial rapid increase in L corresponds to the increase in the average size 
of the polygons. The decrease corresponds to the breakup of polygons into nanoparticles, with 
the long time value representing R. The size distribution of nanoparticles in the final state was 
observed to be monomodal, as shown in Fig. • Also, as seen from Fig. EJB, the value of R 
remains unchanged over long irradiation times and is independent of laser energy E. However, 
increasing E increases the rate at which the patterns evolve, as evident from Fig. |2B and the figure 
in the inset. At the lowest energy of 110 mJ /cm? the time to reach R (n > 3000) was observed 
to be much larger than the 130 m J /cm 2 case ( n < 3000). Since the temporal scale used here is 
coarse, i.e., it is expressed in number of pulses, the precise location of the peak value of L was not 
detectable for every energy value. 

4 Discussion 

It is well accepted that dewetting morphology can progress via three pathways [29]: (i) Homoge- 
neous nucleation and growth, where holes appear spontaneously at random locations and times on 
the surface. Because of the inherent randomness, no characteristic length is present in this type 
of dewetting [20]; (ii) Heterogeneous nucleation and growth due to defects, impurities or other 
experimentally imposed heterogeneities. In this type of dewetting, a characteristic length scale 
could appear at the early stages of dewetting due to ordered nucleation sites. For instance, in ion- 
irradiation induced dewetting, the average molten zone of an ion imposes a characteristic length 
scale in dewetting [23]; and (iii) Thin film hydrodynamic (T.F.H.) instabilities such as the one 
associated with the dewetting of spinodally unstable systems. The resulting patterns are charac- 
terized by a well-defined length scale in the hole spacing and/or size [19]. Since we have ruled 
out spatially ordered heterogeneities on the surface and in the film microstructure, the evidence for 
spatial ordering in Fig. [T] strongly suggests spinodal dewetting. 

We examined whether T.F.H. models developed based on lubrication analysis could help ex- 
plain the time and length scales of the observed dewetting patterns. Towards this end, we adopted 
the T.F.H. equations employed by Vrij [13,30] and Sharma and Ruckenstein [14] to study spin- 
odal dewetting in ultrathin polymeric films. Lubrication theory provides a relationship between 
the mass flux (which can be related to temporal variations in film height) and the effective pres- 
sure gradient, which for a thin film could result from perturbations that produce interface curvature, 
and disjoining pressure due to long-range intermolecular forces such as van der Waals interactions. 
For a thin film resting on a substrate contribution to the disjoining pressure can be expressed as 
A/Qirh' 3 where A is the Hamaker constant and h! is the instantaneous film thickness. This approach 
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leads to the dynamical equation describing the film height as a function of time as: 

dh' / 9 . A* , A 

3 ^ = ~ v • h h v-v 2 ti + — vh') (1) 



dt v V V ^ 

where 77 is the viscosity and 7 is the interfacial energy of the film-vacuum interface and A* = 
A/2n. If infinitesimaHy small perturbations h are imposed on an initially flat film of uniform 
thickness h, i.e. h' = h + h, the growth/decay rate of such perturbations can be assessed by using 
linear stability analysis in which one keeps only terms which are linear in h: 

3r^ = -V. ( 7/1 3 V. V 2 /i + 4^ Vh) . (2) 
at \ h J 

Substitution of the normal mode form h = eexp(ik x x + ik y y + at), where x and y denote 
the orthogonal surface coordinates and a represents the growth or decay rate of perturbations with 
the characteristic wavevector k = (k x , k y ), into eq. EJgives the following dispersion relationship 
(k 2 = k.k): 

3 W = --fh 3 k 4 + 4^ 2 - (3) 
h 

It is evident from eq. |3] that a > for perturbations with wavevectors with k 2 < implying 
that such perturbations will be amplified as time progresses. The fastest growing mode k m can 

1 = 0- This gives k 2 m = 5^. 

A = 2n /k m is given by: 



be obtained by letting || = 0. This gives k 2 m = -^h- Hence, the characteristic wavelength 



The time scale r D = 2n/a associated with the growth of perturbations with wavelength A can be 
evaluated by letting k = k m in eq. EJas: 

967r 3 7?7/i 5 

td = ^ • (5) 

A previous investigation by Bischof et al. using single shot ns pulsed laser melting of a multilayer 
metal structure has shown that is possible to capture snapshots of a dewetting instability provided 
the characteristic time scale of the dewetting process is much larger than the liquid melt time 
r m [31,32]. Based on our analytical model the liquid lifetimes were estimated to be 2 < tl < 12ns 
for the range of experimental parameters reported here [27]. Further, using eq. |5]we estimated the 
dewetting time scale at the melting point T me i t of Co using 7 Co (T m ) = 1.88 J/m 2 , viscosity 
r] Co {T m ) = 4.46 x 1CT 3 Pa — s, and a typical range of the Hamaker constant A for metals from 
IQ-20 < A < 1CT 18 J. We found that the lower limit is m ~ 25ns for a lnm thick film and 
this time increases rapidly with increasing h. Therefore, it is reasonable to expect that pulsed 



6 



laser melting used in this experiment could capture snapshots of spinodal dewetting. We further 
suggest that once the dewetting instability is initiated, continued laser pulsing fosters the instability 
eventually leading to a robust final state whose length scale, as suggested by eq. El should be 
independent of the laser parameters. This effect is observed in Fig. EJB, where despite a change in 
laser E and/or n, an identical NN spacing R is observed in the stable nanoparticle state. 

To provide further quantitative evidence for spinodal dewetting, we have investigated the ex- 
perimental variation of the NN particle spacing R with h. The origin of the space-filling polygonal 
shapes can be explained by the presence of the preferred spinodal wavenumber k ~ \, which is a 
scalar measure of the wavevectors k x and k y in the two orthogonal directions along the undisturbed 
planar film, defined as k = ^Jk 2 + ky. Since there are no external length scales in the x and y di- 
rections, one could have various combinations of k x and k y leading to the observed shapes [33]. 
Therefore the result of the linear spinodal theory (eq. |4j) suggests that the polygon characteristics 
such as size, center-to-center spacings and the vertex separations will be proportional to h 2 . As we 
have shown from Fig. d the observed nanoparticles form predominantly at the polygon vertices. 
Therefore it is clear that if the nanoparticle spacing is a result of spinodal dewetting, then R should 
be proportional to A and also vary as h 2 . In Fig. EfD), the experimentally observed variation of 
R (solid circles) with h is plotted. A power law fit to the variation in R yields R = 25.7/i L98±a3 
where the exponent of h is in excellent agreement with the prediction of eq. E|of 2. The NN 
spacing varied from ~ 50 to 1000 ram for the films studied. From the prefactor of 25.7 ram -1 , 
the Hamaker constant was estimated as A = 14.1 x 10~ 19 J (assuming the proportionality factor 
correlating R to A is 1). This value is of the right sign and order of magnitude for metals in gen- 
eral [34]. We also evaluated the average diameter D of the nanoparticles as a function of h (Fig. 
13) open squares). D varied from ~ 30 to 250 ram for the thickness range studied here. Using 
a volume conservation argument D should vary with h as D ~ h 5 ^ 3 [18]. The power-law fit to 
our experimental data (Fig. |2t>) gave D = 14.1/i L6±0 3 where the exponent is once again in good 
agreement with the theoretical prediction. These results clearly support a dewetting mechanism 
caused by a thin film hydrodynamic instability and also provides evidence that the linear spinodal 
theory is applicable under experimental conditions in which multiple cycles of film melting and 
solidification are present. 

The break-up of the polygonal structures into droplets requires some elaboration. Previous 
studies of break-up of cellular structures into particles have attributed the process to a Rayleigh 
instability [35] as a result of which the particle diameter will follow the trend D ~ h 3 ^ 2 [18,36]. 
However, our interpretation of the experimental observations reported here is that particle forma- 
tion occurs via progression of the spinodal dewetting instability and through capillary effects. The 
primary evidence that supports this interpretation is that the nanoparticles predominantly form at 
the junctions of the polygons, as compared to the classical string-like particle distribution result- 
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ing from the Rayleigh instability of the polygon sides [18]. Capillary effects resulting from the 
difference in curvature of the spheroidal particles at the junctions and the polygon walls will also 
act to drive the liquid towards the junction, a situation that is relatively well understood from the 
classical problem of dragging of a liquid film by a moving plate [37]. 

For the spinodal process it is possible to rationalize why R is invariant to laser energy E despite 
the increase in the rate of dewetting a oc L/n, as observed from Fig. EtB) for early stages (small 
values of n). The increase in rate with increasing E is due to two effects: one is the increase in 
liquid lifetime r L which will prolong the dewetting process within a given number of pulses; the 
second is the increase in liquid temperature T L above T m which will increase the rate as: 

Act(AT) = -a{T m )[-±—% + -±rr%]lST 
7(T m ) dT r)[T m ) dT 

where the Hamaker constant was treated as T-independent based on its extremely weak In- 
dependence [38]. Since the the T-dependent material parameters for Co of dj/dT = —0.5 x 
1CT 3 J/m 2 — s ; and dr//dT = —6.3 x 10~ 6 Pa — s/K both decrease with increasing T, therefore, 
the rate increases with increasing T, as qualitatively observed from the inset of Fig. EtB). This 
increase will be reinforced by that produced by the longer melt lifetime. A similar analysis shows 
that the dewetting length scales L oc A will decrease with increasing AT as: 

AL(AT) = i^I^AT = -2.66 x lCT 4 L(T m )AT 

This decrease in L with increasing T L can be estimated to be ~ —13% for AT = 500 K, which 
is a relatively small change and will be undetectable given the larger spread in the experimental 
values of R of ±20% as measured from the g(k). Therefore, for conditions under which the net 
Co evaporation is small, the changes in R will be experimentally indistinguishable. 



5 Conclusion 

In conclusion, we studied pattern evolution in nanoscopic Co films (1 < h < 8 nm) under irradi- 
ation by multiple 9 ns laser pulses. Laser melting initiates a dewetting instability and subsequent 
irradiation fosters this dewetting to a stable state of spatially ordered nanoparticles. The dewetting 
instability appears to be spinodal in nature based on the observations of a characteristic length 
scale at all stages of dewetting and the ~ h 2 behavior of the NN spacing of the nanoparticles. 
To our knowledge this is the first experimental evidence for: ( i) formation of metal nanoparticles 
with spatial order via spinodal dewetting; ( ii) evidence for linear spinodal theory under multiple 
cycles of phase change; and (Hi) evidence for robustness of the dewetting process with respect to 
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experimental conditions. The short range order and particle size produced by the pattern selection 
mechanisms can be manipulated by employing films with varying thickness. This offers oppor- 
tunities to develop cost-effective and easily controllable nanomanufacturing processes for metal 
nanoarrays. 

RK acknowledges support by the National Science Foundation through grant # DMI-0449258. 
The authors acknowledge invaluable comments provided by Prof.'s Ken Kelton and Linday Greer 
during the preparation of this manuscript. 
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Figure [TJ Scanning electron micrgraphs of dewetting patterns as a function of increasing 
number of laser pulses n for a 3.8 nm thick Co film at a laser energy density of 120 mJ / cm 2 . 
(A) Polygons after n = 10; (B) Coalescing polygons with evidence for particles at the 
junctions after n = 500; (C) Nanoparticles with some remaining evidence of the polygon 
edges after n = 3000; and (D) The final nanoparticle morphology after n = 10500. The 
inset of each figure shows the corresponding power spectra for each pattern. 

Figure I2J Evidence for the characteristic length scales and their behavior for the dewetting 
process. (A) The g(k) (symbols) and fits (lines) for the patterns of Fig. [lj (B) Dependence 
of L on n at various laser energies for the 3.8 nm film. The dashed lines are only intended to 
be guides for the eye. The inset shows the change in L for various energies at early stages, as 
measured after 10 pulses. (C) The monomodal size distribution of nanoparticles in the final 
state for the 3.8 nm film. (D) Nearest-neighbour spacing R (solid circles) and Diameter D 
(open squares) of the nanoparticles vs. film thickness h and the least-squares fit. 
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